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I. INTRODUCTION 



Scalar mesons are still a puzzle in light meson spec- 
troscopy. New claims for the existence of broad states 
close to threshold such as k(800) [l[ and /o(600) Q, have 
reopened discussion about the composition of the ground 
state J PC = ++ nonet, and about the possibility that 
states such as the ao(980) or / (980) may be 4-quark 
states, due to their proximity to the KK threshold Q- 
This hypothesis can be tested only through accurate mea- 
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surements of the branching fractions and the couplings 
to different final states. It is therefore important to have 
precise information on the structure of the 7T7t and KK 
S- waves. In this context, mesons can shed light on 
the structure of the scalar amplitude coupled to ss. The 
7T7T 5-wave has been already extracted from BABAR data 
in a Dalitz plot analysis of Df —> n + ir~TT + 4]. The 
understanding of the KK 5-wave is also of great impor- 
tance for the precise measurement of CP-violation in B s 
oscillations using B s — > J/0 [j| @ . 

This paper focuses on the study of meson decay 
to K + K~tt + Q. Dalitz plot analyses of this decay mode 
have been performed by the E687 and CLEO collabora- 
tions using 700 events [8|, and 14400 events |9| respec- 
tively. The present analysis is performed using about 
100, 000 events. 

The decay — > 07r + is frequently used in particle 
physics as the reference mode for decay. Previous 
measurements of this decay mode did not, however, ac- 
count for the presence of the KK S-w&ve underneath the 
peak. Therefore, as part of the present analysis, we 
obtain a precise measurement of the branching fraction 
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B(Df 4>tt+) relative to B(Df K + K~ir + ). 

Singly Cabibbo-suppressed (SCS) and doubly 
Cabibbo-suppressed (DCS) decays play an important 
role in studies of charmed hadron dynamics. The naive 
expectations for the rates of SCS and DCS decays are of 
the order of tan 2 8c and tan 4 8c, respectively, where Qq 
is the Cabibbo mixing angle. These rates correspond to 
about 5.3% and 0.28% relative to their Cabibbo-favored 
(CF) counterpart. Due to the limited statistics in past 
experiments, branching fraction measurements of DCS 
decays have been affected by large statistical uncertain- 
ties [l(l. A precise measurement of ^7^4—^-^ ~ A nas 

1 — J L B(DJ—>-K+K~w +j 

been recently performed by the Belle experiment |ll| . 
In this paper we study the Df decay 

L>+ -> K+K-JT+ (1) 

and perform a detailed Dalitz plot analysis. We then 
measure the branching ratios of the SCS decay 

Df -> K + K~K + (2) 

and the DCS decay 

D+ -> K+K+tt- (3) 

relative to the CF channel ((T|). The paper is organized 
as follows. Section |TT] briefly describes the BABAR detec- 
tor, while Sec. IIIII gives details of event reconstruction. 
Section [TV] is devoted to the evaluation of the selection 
efficiency. Section [V] describes a partial wave analysis 
of the K + K~ system, the evaluation of the Df — > 4>ir + 
branching fraction and the KK 5-wave parametrization. 
Section I VII deals with the description of the Dalitz plot 
analysis method and background description. Results 
from the Dalitz plot analysis of Df — > K + K~tt + are 
given in Sec. lVIII The measurements of the Df SCS and 
DCS branching fractions are described in Sec. I Villi while 
Sec. IIXI summarizes the results. 



II. THE BABAR DETECTOR AND DATASET 

The data sample used in this analysis corresponds to 
an integrated luminosity of 384 fb _1 recorded with the 
BABAR detector at the SLAC PEP-II collider, operated at 
center-of-mass (cm.) energies near the T(4S) resonance. 
The BABAR detector is described in detail elsewhere [l2j ■ 
The following is a brief summary of the components im- 
portant to this analysis. Charged particle tracks are de- 
tected, and their momenta measured, by a combination 
of a cylindrical drift chamber (DCH) and a silicon vertex 
tracker (SVT), both operating within a 1.5 T solenoidal 
magnetic field. Photon energies are measured with a 
CsI(Tl) electromagnetic calorimeter (EMC). Information 
from a ring-imaging Cherenkov detector (DIRC), and 
specific energy-loss measurements in the SVT and DCH, 
are used to identify charged kaon and pion candidates. 



III. EVENT SELECTION AND D+ K + K~-n + 
RECONSTRUCTION 

Events corresponding to the three-body Df —> 
K + K~ir + decay are reconstructed from the data sam- 
ple having at least three reconstructed charged tracks 
with net charge ± 1. We require that the invariant 
mass of the K + K~-k + system lie within the mass interval 
[1.9-2.05] GeV/c 2 . Particle identification is applied to the 
three tracks, and the presence of two kaons is required. 
The efficiency that a kaon is identified is 90% while the 
rate that a kaon is misidentified as a pion is 2%. The 
three tracks are required to originate from a common ver- 
tex, and the x 2 fit probability (Pi ) must be greater than 
0.1%. We also perform a separate kinematic fit in which 
the Df mass is constrained to its known value [10( . This 
latter fit will be used only in the Dalitz plot analysis. 

In order to help in the discrimination of signal from 
background, an additional fit is performed, constraining 
the three tracks to originate from the e + e~ luminous re- 
gion (beam spot). The x 2 probability of this fit, labeled 
as P 2 , is expected to be large for most of the background 
events, when all tracks originate from the luminous re- 
gion, and small for the Df signal, due to the measurable 
flight distance of the latter. 

The decay 

£> s *(2112)+ -> £>+7 (4) 

is used to select a subset of event candidates in order 
to reduce combinatorial background. The photon is re- 
quired to have released an energy of at least 100 MeV 
into the EMC. We define the variable 

Am = m(K+K-TT + j) - m(K + K~n + ) (5) 

and require it to be within ±2<t d *+ with respect to 
Am D .+ where Am D «+ = 144.94 ± 0.03 stat MeV/c 2 and 
cr D »+ = 5.53±0.04 sta t MeV/c 2 are obtained from a Gaus- 
sian fit of the Am distribution. 

Each Df candidate is characterized by three variables: 
the cm. momentum p* in the e + e~ rest frame, the differ- 
ence in probability P1 — P2, and the signed decay distance 
dxv = t Pa 7 where d is the vector joining the beam spot 

to the Df decay vertex and p xy is the projection of the 
Df momentum on the xy plane. These three variables 
are used to discriminate signal from background events: 
in fact signal events are expected to be characterized by 
larger values of p* [r|, due to the jet-like shape of the 
e + e~ — » cc events, and larger values of d xy and P\ — P2, 
due to the measurable flight distance of the Df meson. 

The distributions of these three variables for signal 
and background events are determined from data and 
are shown in Fig. [TJ The background distributions are 
estimated from events in the Df mass-sidebands, while 
those for the signal region are estimated from the 
signal region with sideband subtraction. The normalized 
probability distribution functions (PDFs) are then com- 
bined in a likelihood-ratio test. A selection is performed 
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3 4 5 -0.1 0.1 0.2 0.3 0.5 1 



p* (GeV/c) d xy (cm) P, - P 2 

FIG. 1: Normalized probability distribution functions for signal (solid) and background events (hatched) used in a likelihood- 
ratio test for the event selection of Df — > K + K~ n + : (a) the center of mass momentum p* , (b) the signed decay distance d xy 
and (c) the difference in probability Pi — P-2- 



on this variable such that signal to background ratio is 
maximized. Lower sideband, signal and upper sideband 
regions are denned between [1.911 - 1.934] GeV/c 2 , [1.957 
- 1.980] GeV/c 2 and [2.003 - 2.026] GeV/c 2 , respectively, 
corresponding to (— 10er, — 6a), (— 2a, 2a) and (6a, 10<r) 
regions, where a is estimated from the fit of a Gaussian 
function to the Df lineshape. 

We have examined a number of possible background 
sources. A small peak due to the decay D* + — > n + D° 
where D° — > K + K~ is observed. A Gaussian fit to this 
K + K~ spectrum gives a D o^, K + K - — 5.4 MeV/c 2 . For 
events within 3.5aD°->K+K- of the D° mass, we plot 
the mass difference Am(K + K~ir + ) — m(K + K~tt + ) — 
m(K + K~) and observe a clean D* + signal. We remove 
events that satisfy Am(K+K-TT+) < 0.15 GeV/c 2 . The 
surviving events still show a D° — > K + K~ signal which 
does not come from this D* + decay. We remove events 
that satisfy m(K + K~) >1.85 GeV/c 2 . 

Particle misidentification, in which a pion 7Lt. is 
wrongly identified as a kaon, is tested by assigning the 
pion mass to the K + . In this way we identify the back- 
ground due to the decay D + — > K~-k + ii + which, for the 
most part, populates the higher mass Df — > K + K~n + 
sideband. However, this cannot be removed without bi- 
asing the Dalitz plot, and so this background is taken 
into account in the Dalitz plot analysis. 

We also observe a clean peak in the distribution of the 
mass difference 77j(.ftT - 7r^ lis 7r + ) — m{K~ it^-^) . Combin- 
ing m(.ftT - 7r^ is ) with each of the n° meson candidates 
in the event, we identify this contamination as due to 
D*+ -> tt+D"(-^ K-jt+tt ) with a missing tt°. We re- 
move events that satisfy m{K~-K^ nis -n + ) — m{K~ 7r+ is ) < 
0.15 GeV/c 2 . Finally, we remove the Df candidates 
that share one or two daughters with another Df candi- 
date; this reduces the number of candidates by 1.8%, 
corresponding to 0.9% of events. We allow there to 



be two or more non-overlapping multiple candidates in 
the same event. The resulting K + K~ir + mass dis- 
tribution is shown in Fig. Ufa). This distribution is 
fitted with a double-Gaussian function for the signal, 
and a linear background. The fit gives a mass of 
1968.70 ±0.02 stat MeV/c 2 , err = 4.96±0.06 sta t MeV/c 2 , 
02/01 = 1.91 ± 0.06 s tat where a\ (a 2) is the standard 
deviation of the first (second) Gaussian, and errors are 
statistical only. The fractions of the two Gaussians are 
f ai = 0.80 ± 0.02 and f a2 = 0.20 ± 0.02. The signal 
region is defined to be within ±2a D + of the fitted mass 

value, where a D + = \/ f^a'f + f a ^a\ = 6.1 MeV/c 2 is 
the observed mass resolution (the simulated mass resolu- 
tion is 6 MeV/c 2 ) . The number of signal events in this 
region (Signal), and the corresponding purity (defined as 
Signal/(Signal+Background)), are given in Table Q] 

TABLE I: Yields and purities for the different Df decay 
modes. Quoted uncertainties are statistical only. 



D^t decay mode 


Signal yield 


Purity (%) 


K+K--R+ 


96307 ± 369 


95 


K+K-K+ 


748 ± 60 


28 


K+K+ty- 


356 ± 52 


23 



For events in the Df — > K + K~tt + signal region, we 
obtain the Dalitz plot shown in Fig. &h). For this dis- 
tribution, and for the Dalitz plot analysis (Sec. IVI[) . we 
use the track parameters obtained from the Df mass- 
constrained fit, since this yields a unique Dalitz plot 
boundary. 

In the K + K~ threshold region, a strong 0(1020) signal 
is observed, together with a rather broad structure. The 
/o(980) and cio(980) 5-wave resonances are, in fact, close 
to K + K~ threshold, and might be expected to contribute 
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FIG. 2: (a) K + K tt + mass distribution for the D + analysis sample; the signal region is as indicated, (b) D + — > K + K 7r H 
Dalit z plot. 



in the vicinity of the 0(1020). A strong K*(892)° signal 
can also be seen in the K~n + system, but there is no 
evidence of structure in the K + ir + mass. 



IV. EFFICIENCY 

The selection efficiency for each D + decay mode ana- 
lyzed is determined from a sample of Monte Carlo (MC) 
events in which the D + decay is generated according to 
phase space (i.e. such that the Dalitz plot is uniformly 
populated). The generated events are passed throu gh a 
detector simulation based on the Geant4 toolkit (14| . 
and subjected to the same reconstruction and event se- 
lection procedure as that applied to the data. The dis- 
tribution of the selected events in each Dalitz plot is 
then used to determine the reconstruction efficiency. The 
MC samples used to compute these efficiencies consist 
of 4.2 xlO 6 generated events for Df — > K + K~tt + and 
Df ->■ K+K+tt-, and 0.7 xlO 6 for Df K+R-R+ . 

For D + — > K + K~ir + , the efficiency distribution is fit- 
ted to a third-order polynomial in two dimensions using 
the expression: 

r](x, y) = ao + a\x' + a^x' 2 + a^y' 2 + a^x'y' 

+a 6 x' 3 + a 7 y' 3 (6) 

where x — m 2 (K + K~), y — m 2 (K~ir + ), x 1 — x—2, and 
y' = y — 1.25. Coefficients consistent with zero have been 
omitted. We obtain a good description of the efficiency 
with X 2 /NDF = 1133/(1147 - 7) = 0.994 (NDF = 
Number of Degrees of Freedom) . The efficiency is found 
to be almost uniform in K~n + and K + K~ mass, with 
an average value of ~ 3.3% (Fig. [3]). 



V. PARTIAL WAVE ANALYSIS OF THE K+ K~ 
AND K~-k + THRESHOLD REGIONS 

In the K+K- threshold region both ao(980) and 
io(980) can be present, and both resonances have very 
similar parameters which suffer from large uncertainties. 
In this section we obtain model-independent information 
on the K + K~ 5-wave by performing a partial wave anal- 
ysis m the K + K~ threshold region. 

Let N be the number of events for a given mass interval 
/ = [niK+K- 1 m K+K~ + drriK+K-]- We write the corre- 
sponding angular distribution in terms of the appropriate 
spherical harmonic functions as 

k=0 

where L = 2£ max , and t max is the maximum orbital an- 
gular momentum quantum number required to describe 
the K + K~ system at m K + K - (e.g. £ max = 1 for an 
S-, "P-wave description); 6 is the angle between the K + 
direction in the K + K~ rest frame and the prior direc- 
tion of the K + K~ system in the D + rest frame. The 
normalizations are such that 

/l r 
y fe °(cos6')Y 7 (cos6i)dcos6' = (8) 
.i 2tt 

and it is assumed that the distribution , „ has been 

d cos & 

efficiency-corrected and background-subtracted. 

Using this orthogonality condition, the coefficients in 
the expansion are obtained from: 

(Y k 0) = lY^^cose (9) 



9 



■ 0.05 




m 2 (K + K") (GeV 2 /c 4 ) 



m 2 (K + K") (GeV 2 /c 4 ) 



m 2 (K'7i + ) (GeV7c ) 

FIG. 3: (a) Dalitz plot efficiency map; the projection onto (b) the K + K~ and (c) the K~tt + axis. 
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FIG. 4: K + K mass spectrum in the threshold region weighted by (a) Yo°, (b) Y®, and (c) Y®, corrected for efficiency and 
phase space, and background-subtracted. 



where the integral is given, to a good approximation, by 
J2n=i Yk(cos9 n ), where 9 n is the value of 9 for the n-th 
event. 

Figure 2] shows the K + K~ mass spectrum 
up to 1.5GeV/c 2 weighted by Y£(cos9) = 
y / {2k + 1)/4tt P k {cos9) for k = 0,1 and 2, where 
Pfe is the Legendre polynomial function of order k. 
These distributions are corrected for efficiency and 
phase space, and background is subtracted using the 
sidebands. 

The number of events N for the mass interval I can be 
expressed also in terms of the partial-wave amplitudes de- 
scribing the K + K~ system. Assuming that only S- and 
P-wave amplitudes are necessary in this limited region, 
we can write: 



= 2n\SY°(cos9)+VY°(cos9)\ 2 . (10) 



By comparing Eq. ([7]) and Eq. (JTUJ) !l5[, we obtain: 

V^(Y°) = |S| 2 + |P| 2 

V^(Y°) = 2|«S||P|cos0 S7 > (11) 

V5F<1?) = ±\V\ 2 

where §sv — <Ps ~ 4>v is the phase difference be- 
tween the S- and P-wave amplitudes. These equa- 
tions relate the interference between the S- wave (/o(980), 
and/or ao(980), and/or nonresonant) and the P-wave 
(0(1020)) to the prominent structure in (Y?) (Fig.QJb)). 
The (Yi°) distribut ion shows the same behavior as for 
D+ ->• K + K-e + v e decay [Hj]. The (F 2 °) distribution 
(Fig. IHc)), on the other hand, is consistent with the 
0(1020) lineshape. 

The above system of equations can be solved in each in- 
terval of K + K~ invariant mass for \S\, |P|, and 4>sV: and 
the resulting distributions are shown in Fig. [SJ We ob- 
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serve a threshold enhancement in the S- wave (Fig. 03a)), 
and the expected 0(1020) Breit-Wigner (BW) in the V- 
wave (Fig. (5][b)). We also observe the expected S-V rel- 
ative phase motion in the 0(1020) region (Fig. G3c)). 



A. "P-wave/<S-wave ratio in the 0(1020) region 

The decay mode Df — ► 0(lO2O)7r + is used often as the 
normalizing mode for decay branching fractions, typ- 
ically by selecting a K + K~ invariant mass region around 
the 0(1020) peak. The observation of a significant iS- 
wave contribution in the threshold region means that this 
contribution must be taken into account in such a proce- 
dure. 

In this section we estimate the P-wave/ 5- wave ratio 
in an almost model-independent way. In fact integrat- 
ing the distributions of y/Air pq' (Yq) and \/5npq' (Y^) 
(Fig. S| in a region around the 0(1020) peak yields 
f{\S\ 2 + \V\ 2 )pq'dm K+K - and f \V\ 2 pq'dm K+K - respec- 
tively, where p is the K + momentum in the K + K~ rest 
frame, and q' is the momentum of the bachelor n + in the 
rest frame. 

The S-V interference contribution integrates to zero, 
and we define the P-wave and 5-wave fractions as 



f'P— wave 
fs— wave 



J \V\ 2 pq'dm K+K - 
f(\S\* + \V\ 2 )pq'dm K+K - 
J \S\ 2 pq'dm K+K - 



${\S\ 2 + 
1 - fr-y 



\V\ 2 )pq'dm 



k+k- 



(12) 



(13) 



The experimental mass resolution is estimated by com- 
paring generated and reconstructed MC events, and is 
~ 0.5 MeV/c 2 at the mass peak. Table HTl gives the re- 
sulting 5-wave and P-wave fractions computed for three 
K + K~ mass regions. The last column of Table HTl shows 
the measurements of the relative overall rate de- 
fined as the number of events in the K + K~ mass interval 
over the number of events in the entire Dalitz plot after 
efficiency-correction and background-subtraction. 



TABLE II: S-wave and P-wave fractions computed in three 
K + K~ mass ranges around the 0(1020) peak. Errors are 
statistical only. 



m K+K - (MeV/c^) f s 



(%) fv- 



m 



1019.456 ±5 3.5 ± 1.0 96.5 ± 1.0 29.4 ± 0.2 
1019.456 ± 10 5.6 ± 0.9 94.4 ± 0.9 35.1 ± 0.2 
1019.456 ± 15 7.9 ± 0.9 92.1 ± 0.9 37.8 ± 0.2 



B. <S-wave parametrization at the K + K threshold 

In this section we extract a phenomenological descrip- 
tion of the 5-wave assuming that it is dominated by the 
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FIG. 5: Squared (a) 5- and (b) P-wave amplitudes; (c) the 
phase difference (f>sv, (d) 4>s obtained as explained in the 
text. The curves result from the fit described in the text. 
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/o(980) resonance while the P-wave is described entirely 
by the 0(1020) resonance. We also assume that no other 
contribution is present in this limited region of the Dalitz 
plot. We therefore perform a simultaneous fit of the three 
distributions shown in Figs. Eta), (b), and (c) using the 
following model: 



dN 



5 2 



dm K+K - 

dN V 2 
dm K+K - 



I C/o (980)^/o (980) I 2 



= \C A A, 



(14) 



: arg(A fo{980) e' lS ) - arg{A^) 



dm K+K - 

where Cj,, Cy (g S0 ), and 5 are free parameters and 

F r Fn 



Aa, = 



UK: 



to 2 — im^T 



x Apq 



(15) 



is the spin 1 relativistic BW parametrizing the 0(1020) 
with r expressed as: 



r = r r 



2,7+1 



M r 
m 



F 



(16) 



Here q is the momentum of the bachelor 7r + in the 
K + K~ rest frame. The parameters in Eqs. (fT5j) and 
(fTB|) are defined in Sec. IVII below. 



For A 



/o(980) 



we first tried a coupled channel BW 



(Flatte) amplitude 17]. However we find that this 
parametrization is insensitive to the coupling to the 
7T7T channel. Therefore we empirically parametrize the 
/o(980) with the following function: 



.4 



1 



/o(980) 



(17) 



where Pkk 
values: 



- to 2 - im T p KK 
2p/m, and obtain the following parameter 



too = (0.922 ± 0.003 sta t) GeV/c 2 
T = (0.24 ±0.08 

slat 

) GeV 



(18) 



The errors are statistical only. The fit results are super- 
imposed on the data in Fig. [5] 

In Fig. [SJc), the S-V phase difference is plotted twice 
because of the sign ambiguity associated with the value 
of 4>sv extracted from cos 4>sv- We can extract the mass- 
dependent /o(980) phase by adding the mass-dependent 
0(1020) BW phase to the 4>sv distributions of Fig.[5jc). 
Since the K + K~ mass region is significantly above the 
/o(980) central mass value of Eq. (fT8|) . we expect that 
the 5-wave phase will be moving much more slowly in 
this region than in the 0(1020) region. Consequently, 
we resolve the phase ambiguity of Fig. [SJc) by choosing 
as the physical solution the one which decreases rapidly 
in the 0(1020) peak region, since this reflects the rapid 
forward BW phase motion associated with a narrow res- 
onance. The result is shown in Fig. [HJd), where we see 



TABLE III: 5- and P-wave squared amplitudes (in arbitrary 
units) and 5-wave phase. The 5-wave phase values, corre- 
sponding to the mass 0.988 and 1.116 GeV/c 2 , are missing 
because the (Y®) distribution (Fig. [He)) goes negative or 
\cos 4>sv\ > 1 and so Eqs. ()lip cannot be solved. Quoted 
uncertainties are statistical only. 



m K+K - 
(GeV/c 2 ) 



(arbitrary units) 



w 

(arbitrary 



units) 



4>s 

(degrees) 



0.988 


22178 


± 


3120 


-133 


± 


2283 








0.992 


18760 


± 


1610 


2761 


± 


1313 


92 


± 


5 


0.996 


16664 


± 


1264 


1043 


± 


971 


84 


± 


7 


1 


12901 


± 


1058 


3209 


± 


882 


81 


± 


4 


1.004 


13002 


± 


1029 


5901 


± 


915 


82 


± 


3 


1.008 


9300 


± 


964 


13484 


± 


1020 


76 


± 


3 


1.012 


9287 


± 


1117 


31615 


± 


1327 


80 


± 


2 


1.016 


6829 


± 


1930 


157412 


± 


2648 


75 


± 


8 


1.02 


11987 


± 


2734 


346890 


± 


3794 


55 


± 


6 


1.024 


5510 


± 


1513 


104892 


± 


2055 


86 


± 


5 


1.028 


7565 


± 


952 


32239 


± 


1173 


75 


± 


2 


1.032 


7596 


± 


768 


15899 


± 


861 


74 


± 


2 


1.036 


6497 


± 


658 


10399 


± 


707 


77 


± 


2 


1.04 


5268 


± 


574 


7638 


± 


609 


72 


± 


3 


1.044 


5467 


± 


540 


5474 


± 


540 


72 


± 


3 


1.048 


5412 


± 


506 


4026 


± 


483 


72 


± 


3 


1.052 


5648 


± 


472 


2347 


± 


423 


71 


± 


3 


1.056 


4288 


± 


442 


3056 


± 


421 


70 


± 


3 


1.06 


4548 


± 


429 


1992 


± 


384 


73 


± 


3 


1.064 


4755 


± 


425 


1673 


± 


374 


70 


± 


4 


1.068 


4508 


± 


393 


1074 


± 


334 


75 


± 


4 


1.072 


3619 


± 


373 


1805 


± 


345 


75 


± 


4 


1.076 


4189 


± 


368 


840 


± 


312 


70 


± 


5 


1.08 


4215 


± 


367 


770 


± 


297 


71 


± 


5 


1.084 


3508 


± 


345 


866 


± 


294 


71 


± 


5 


1.088 


3026 


± 


322 


929 


± 


285 


75 


± 


4 


1.092 


3456 


± 


309 


79 


± 


240 


37 


± 


90 


1.096 


2903 


± 


300 


488 


± 


256 


75 


± 


6 


1.1 


2335 


± 


282 


885 


± 


248 


68 


± 


5 


1.104 


2761 


± 


284 


341 


± 


231 


57 


± 


10 


1.108 


2293 


± 


273 


602 


± 


231 


77 


± 


5 


1.112 


1913 


± 


238 


269 


± 


186 


74 


± 


8 


1.116 


2325 


± 


252 


57 


± 


198 








1.12 


1596 


± 


228 


308 


± 


194 


78 


± 


7 


1.124 


1707 


± 


224 


233 


± 


188 


67 


± 


10 


1.128 


1292 


± 


207 


270 


± 


176 


66 


± 


9 


1.132 


969 


± 


197 


586 


± 


172 


60 


± 


6 


1.136 


1092 


± 


196 


553 


± 


170 


67 


± 


6 


1.14 


1180 


± 


193 


316 


± 


167 


48 


± 


11 


1.144 


1107 


± 


187 


354 


± 


170 


68 


± 


8 


1.148 


818 


± 


178 


521 


± 


164 


64 


± 


7 



that the 5-wave phase is roughly constant, as would be 
expected for the tail of a resonance. The slight decrease 
observed with increasing mass might be due to higher 
mass contributions to the 5-wave amplitude. The values 
of |5| 2 (arbitrary units) and phase values are reported in 
Table IhII together with the corresponding values of 1 7^ | 2 - 
In Fig. [6ja) we compare the 5-wave profile from this 
analysis with the 5-wave intensity values extracted from 
Dalitz plot analyses of L»° K 6 K+R- and D° 
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FIG. 6: (a) Comparison between KK 5-wave intensities from different charmed meson Dalitz plot analyses, (b) Comparison 
of the KK 5-wave intensity from Df — > K + K~ir + with the tt + -k~ 5-wave intensity from Df — > 7r + 7r~7r + . 



|19j. The four distributions are normalized in 
the region from threshold up to 1.05 GeV/c 2 . We ob- 
serve substantial agreement. As the ao(980) and /o(980) 
mesons couple mainly to the uu/dd and ss systems re- 
spectively, the former is favoured in D° —> K°K + K~ 
and the latter in Df —> K + K~tt + . Both resonances 
can contribute in D° — > K + K~ir°. We conclude that 
the 5-wave projections in the KK system for both res- 
onances are consistent in shape. It has been suggested 
that this feature supports the h ypo thesis that the ao (980) 
and /o(980) are 4-quark states [20j . We also compare the 
5-wave profile from this analysis with the tt + tt~ 5-wave 
profile extracted from BABAR data in a Dalitz plot anal- 
ysis of D+ 7r+7r-7r+ [4] (Fig. E^b)). The observed 
agreement supports the argument that only the /o(980) 
is present in this limited mass region. 



C. Study of the K 7r + S-wave at threshold 

We perform a model-independent analysis, similar to 
that described in the previous sections, to extract the 
Ktt 5-wave behavior as a function of mass in the thresh- 
old region up to l.lGeV/c 2 . Figure [7] shows the K~n + 
mass spectrum in this region, weighted by Y®(cos9) = 
y/ {2k + l)/47rP fe (cos#), with k = 0, 1 and 2, corrected 
for efficiency, phase space, and with background from 
the Df sidebands subtracted; 8 is the angle between the 
K~ direction in the K~n + rest frame and the prior di- 
rection of the K~tt + system in the rest frame. We 
observe that (Y§) an d (^2°) show strong K* (892)° reso- 
nance signals, and that the (Y^'S moment shows evidence 
for S-V interference. 



We use Eqs. dTTJ) to solve for |5| and \V\. The re- 
sult for the 5-wave is shown in Fig. [TJd). We observe 
a small 5-wave contribution which does not allow us to 
measure the expected phase motion relative to that of 
the K* (892)° resonance. Indeed, the fact that |5| 2 goes 
negative indicates that a model including only 5- and V- 
wave components is not sufficient to describe the i^ _ 7r + 
system. 



VI. DALITZ PLOT FORMALISM 

An unbinned maximum likelihood fit is performed in 
which the distribution of events in the Dalitz plot is used 
to determine the relative amplitudes and phases of inter- 
mediate resonant and nonresonant states. 

The likelihood function is written as: 



n=l 



Us ■ v(x, y) 
(i - /«*) 



Y^ij c i c*A i (x,y)A*(x,y) 



E 4 hlBi 
where: 

• N is the number of events in the signal region; 



(19) 



• x 



m 



'(K+K-) and y = m 2 (K-n+) 



• /sig is the fraction of signal as a function of the 
K + K ~7T + invariant mass, obtained from the fit to 
the K + K~ir + mass spectrum (Fig. [2a)); 

• r)(x, y) is the efficiency, parametrized by a 3 rd order 
polynomial fSec. IIVI) : 
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• the Ai(x, y) describe the complex signal amplitude 
contributions; 

• the Bi(x,y) describe the background probability 
density function contributions; 

• ki is the magnitude of the i-th component for the 
background. The ki parameters are obtained by 
fitting the sideband regions; 

• = J A i {x,y)A*{x 1 y)r](x 1 y)dxdy and 
Ig i = J Bi(x, y)dxdy are normalization integrals. 
Numerical integration is performed by means of 
Gaussian quadrature |21| ; 

• Ci is the complex amplitude of the «-th component 
for the signal. The Cj parameters are allowed to 
vary during the fit process. 

The phase of each amplitude (i.e. the phase of the 
corresponding Cj) is measured with respect to the 
K + K* (892)° amplitude. Following the method described 
in Ref. [22}, each amplitude Ai{x, y) is represented by the 
product of a complex BW and a real angular term T de- 
pending on the solid angle Q: 



A(x,y) = BW(m) x T{Q). 



(20) 



[M 



AB 



M\ 



Mlf-AM\Ml 



PAB = 2M^ ' (23) 

The form factors F r and Fjj attempt to model the un- 
derlying quark structure of the parent particle and the in- 
termediate resonances. We use the Blatt-Weisskopf pen- 
etration factors [23| (Table HV)) , that depend on a single 
parameter R representing the meson "radius". We as- 
sume R D + = 3GeV _1 for the D s and R r = 1.5 GeV" 1 
for the intermediate resonances; qAB is the momentum 
of the bachelor C in the AB rest frame: 



'Ml 



M\ B) 



V V, LI, • O ALSJ LJ S U 

qAB ^ 1 2M^ ■ ™ 

p r and q r are the values of pab and qAB when jtiab = w r . 



TABLE IV: Summary of the Blatt-Weisskopf penetration 
form factors. q r and p r are the momenta of the decay particles 
in the parent rest frame. 



Spin 



F r 



F D 



yJl + (RrP r ) 2 

^/l + (R rPA B) 2 



via an intermediate resonance r (D s — > rC, r 
BW(Mab) is written as a relativistic BW: 



AB )l r, y/9+3(RrPr) 2 HRrPr) i V 


/9+3(fl £j + 9r ) 2 +(B i3 + qr ) 4 


^/g+'MRrPAB^ + iRrPAB) 4 ^ 


lg+3(R D +q A B) 2 + (R o+lAB) 4 - 



BW(Mab) 



F t Fd 



M AB 



ir AB M r 



(21) 



where Tab is a function of the invariant mass of system 
AB (Mab), the momentum pab of either daughter in 
the AB rest frame, the spin J of the resonance and the 
mass M r and the width T r of the resonance. The explicit 
expression is: 



Tab = T, 



(if) 



2,7+1 



(Mr_\ p 2 

\Mab) r 



(22) 



The angular terms T{D) are described by the following 
expressions: 

Spin : T(O) =1 



Spin 1 : T(J?) =M% C - M 



AC 



(M 2 



D, 



M*){M%-M\) 



(25) 



Spin 2 : T{Q) =a\ — -0203 
o 
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where: 



is performed using the function: 



ai =M BC - M AC + s 



a 2 =M\ B - 2M 2 Ds 



2M% 



(26) 



«3 =M\ B - 2M\ - 2M B 



M\ B 



Resonances are included in sequence, starting from 
those immediately visible in the Dalitz plot projections. 
All allowed resonances from Ref. [lOj have been tried, 
and we reject those with amplitudes consistent with zero. 
The goodness of fit is tested by an adaptive binning y 2 . 

The efficiency-corrected fractional contribution due to 
the resonant or nonresonant contribution i is defined as 
follows: 



fi = 



|cj| 2 /|^(a;,y)| 2 da;dy 

f\T,j c j A j( x iy)\ 2dxd y' 



(27) 



The fi do not necessarily add to 1 because of interfer- 
ence effects. We also define the interference fit fraction 
between the resonant or nonresonant contributions k and 
I as: 

hi = 2 J^[c k c*A k (x,y)A;(x,y)}dxdy 

I\J2j c j A j( x ^y)\ 2dxd y 

Note that fkk = 2fk- The error on each fi and fki is eval- 
uated by propagating the full covariance matrix obtained 
from the fit. 



A. Background parametrization 

To parametrize the D+ background, we use the 
sideband regions. An unbinned maximum likelihood fit 



JV B 



71=1 



E 2 hiBi 



(29) 



where Nb is the number of sideband events, the ki pa- 
rameters are real coefficients floated in the fit, and the 
Bi parameters represent Breit-Wigner functions that are 
summed incoherently. 

The Dalitz plot for the two sidebands shows the pres- 
ence of 0(1020) and J?*(892)° (Fig.©. There are further 
structures not clearly associated with known resonances 
and due to reflections of other final states. Since they do 
not have definite spin, we parametrize the background 
using an incoherent sum of 5-wave Breit-Wigner shapes. 



VII. DALITZ PLOT ANALYSIS OF 



K+K 



Using the method described in Sec. I VII we perform 
an unbinned maximum likelihood fit to the Dj — > 
K + K~ir + decay channel. The fit is performed in steps, 
by adding resonances one after the other. Most of the 
masses and widths of the resonances are taken from 
Ref. [13 . For the / (980) we use the phenomcnological 
model described in Sec. lVBI The K*(892)° amplitude is 
chosen as the reference amplitude. 

The decay fractions, amplitudes, and relative phase 
values for the best fit obtained, are summarized in Ta- 
bic [V] where the first error is statistical, and the second 
is systematic. The interference fractions are quoted in 
Table I VII where the error is statistical only. We observe 
the following features. 

• The decay is dominated by the K*(892)°K+ and 
0(lO2O)7r+ amplitudes. 

• The fit quality is substantially improved by leaving 
the K*(892)° parameters free in the fit. The fitted 
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TABLE V: Results from the — > K + K ir + Dalitz plot analysis. The table gives fit fractions, amplitudes and phases from 
the best fit. Quoted uncertainties are statistical and systematic, respectively. 



Decay mode Decay fraction 



Amplitude 



Phase (radians) 



K*(892)°K+ 

0(1020) 7T + 

M980) 7T+ 
^(1430)°^ 

/ (1710) 7T + 
/q(1370) 7T + 



47.9 ±0.5 ±0.5 
41.4±0.8±0.5 
16.4±0.7±2.0 
2.4±0.3± 1.0 
1.1±0.1±0.1 
1.1±0.1±0.2 



1. (Fixed) 
1.15±0.01±0.26 
2.67 ±0.05 ±0.20 
1.14 ±0.06 ±0.36 
0.65 ±0.02 ±0.06 
0.46 ±0.03 ±0.09 



0. (Fixed) 
2.89 ±0.02 ±0.04 
1.56 ±0.02 ±0.09 
2.55 ±0.05 ±0.22 
1.36 ±0.05 ±0.20 
-0.45 ±0.11 ±0.52 



Sum 

X 2 /NDF 



110.2 ±0.6 ±2.0 
2843/(2305 - 14) 



= 1.24 



parameters are: 



r 



1 K* (892)° 
if* (892)0 



(895.6 ± 0.2 sta t ± 0.3 sys ) MeV/c 2 
(45.1±0.4 sta t±0.4 sys )MeV 



(30) 



We notice that the width is about 3 MeV lower 
than that in Ref . [10| . However this measurement is 
consistent with results from other Dalitz plot anal- 
yses |{| . 

• The Jo (1370) contribution is also left free in the fit, 
and we obtain the following parameter values: 



m /o( i370) = (1-22 ± 0.01 s tat ± 0.04 sys ) GeV/c 2 
r /o( i370) = (0.21 ± 0.01 s tat ± 0.03 sys ) GeV 



(31) 



These values are within the broad range of values 
measured by other experiments 10]. 

• A nonresonant contribution, represented by a con- 
stant complex amplitude, was included in the fit 
function. However this contribution was found to 
be consistent with zero, and therefore is excluded 
from the final fit function. 

• In a similar way contributions from the if* (1410), 
/o(1500), / 2 (1270), and $(1525) are found to be 
consistent with zero. 

• The replacement of the K*{U30) by the LASS 
parametrization [24[ of the entire Ktt S- wave does 
not improve the fit quality. 

• The fit does not require any contribution from the 
k(800) % 

The results of the best fit ( X 2 /NDF = 2843/(2305 - 
14) = 1.24) are superimposed on the Dalitz plot projec- 
tions in Fig. [§] Other recent high statistics charm Dalitz 
plot analyses at BaBAR [25[ have shown that a significant 
contribution to the x 2 /NDF can arise from imperfec- 
tions in modelling experimental effects. The normalized 
fit residuals shown under each distribution (Fig. [5]) are 
given by Pull = ( AT d ata - N Rt ) / V^data ■ The data are well 



reproduced in all the projections. We observe some dis- 
agreement in the K~ir + projection below 0.5 GeV 2 /c 4 . It 
may be due to a poor parametrization of the background 
in this limited mass region. A systematic uncertainty 
takes such effects in account (Sec. I VII A[) . The missing 
of a Kit iS-wave amplitude in the K~n + low mass region 
may be also the source of such disagreement. 

Another way to test the fit quality is to project the 
fit results onto the {Y£} moments, shown in Fig. 1101 for 
the K + K~ system and Fig. [IT] for the K~ir + system. 
We observe that the fit results reproduce the data pro- 
jections for moments up to k — 7, indicating that the fit 
describes the details of the Dalitz plot structure very well. 
The K ~tt + (i^ ) and (i^ ) moments show activity in the 
K* (892)° region which the Dalitz plot analysis relates to 
interference between the K*(892)°K+ and /o(1710)tt + 
decay amplitudes. This seems to be a reasonable expla- 
nation for the failure of the model- independent K~tt + 
analysis (Sec. IV C[) . although the fit still does not pro- 
vide a good description of the (Y®} and (Y®\ moments 
in this mass region. 

We check the consistency of the Dalitz plot results and 
those of the analysis described in Sec. IV Bl We compute 
the amplitude and phase of the fo (980) /5-wave relative 
to the </>(1020)/'P-wave and find good agreement. 



A. Systematic errors 

Systematic errors given in Table Ivl and in other quoted 
results take into account: 

• Variation of the R r and R D + constants in the Blatt- 
Weisskopf penetration factors within the range [0- 
3] GeV" 1 and [1-5] GeV" 1 , respectively. 

• Variation of fixed resonance masses and widths 
within the ±l<r error range quoted in Ref. [Toj . 

• Variation of the efficiency parameters within ±ltr 
uncertainty. 

• Variation of the purity parameters within ±lcr un- 
certainty. 
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FIG. 9: — >■ K + K 7r + : Dalitz plot projections from the best fit. The data are represented by points with error bars, the fit 
results by the histograms. 



TABLE VI: Fit fractions matrix of the best fit. The diagonal elements /; correspond to the decay fractions in Table FVl The off- 
diagonal elements give the fit fractions of the interference fki ■ The null values originate from the fact that any S-V interference 
contribution integrates to zero. Quoted uncertainties are statitistical only. 



hi (%) 


K*{892)°K+ 


0(1020) 7T+ 


/o (980) 


7T+ 


if 5 (1430)° 


,f (1710) 7T + 


/0 (1370) 7T + 


K*(892)°K+ 


47.9 ± 0.5 


-4.36 ± 0.03 


-2.4 ± 


0.2 


0. 


-0.06 ± 0.03 


0.08 ± 0.08 


0(1020) 7T + 




41.4 ± 0.8 


0. 




-0.7 ± 0.2 


0. 


0. 


/0(980) 7T + 






16.4 ± 


0.7 


4.1 ± 0.6 


-3.1 ± 0.2 


-4.5 ± 0.3 












2.4 ± 0.3 


0.48 ± 0.08 


-0.7 ± 0.1 


/ (1710) Tt + 












1.1 ± 0.1 


0.86 ± 0.06 


/(,(1370) 7T + 














1.1 ± 0.1 
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FIG. 10: K + K~ mass dependence of the spherical harmonic moments, (Y®\, obtained from the fit to the — > K + K~-k + 
Dalitz plot compared to the data moments. The data are represented by points with error bars, the fit results by the histograms. 
The insets show an expanded view of the <£(1020) region. 
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FIG. 11: K 7r + mass dependence of the spherical harmonic moments, (Y°\ obtained from the fit to the Df — > K + K ty + 
Dalitz plot compared to the data moments. The data are represented by points with error bars, the fit results by the histograms. 
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TABLE VII: 
rations. 



Comparison of the fitted decay fractions with the Dalitz plot analyses performed by E687 and CLEO-c collabo- 



Decay mode 



K*{892)°K+ 

0(1020) 7T + 

M980) 7T+ 

/ (1710) 7T+ 
/ (1370) 7T+ 



BABAR 



Decay fraction 
E687 



CLEO-c 



47.9 ± 0.5 ± 0.5 47.8 ± 4.6 ± 4.0 47.4 ± 1.5 ± 0.4 

41.4±0.8±0.5 39.6±3.3±4.7 42.2±1.6±0.3 

16.4±0.7±2.0 11.0±3.5±2.6 28.2±1.9±1.8 

2.4±0.3±1.0 9.3±3.2±3.2 3.9±0.5±0.5 

1.1 ±0.1 ±0.1 3.4±2.3±3.5 3.4±0.5±0.3 

1.1±0.1±0.2 4.3±0.6±0.5 



Sum 

X 2 /NDF 
Events 



110.2±0.6±2.0 

2843 _ i o 
(2305-14) 

96307 ± 369 



111.1 

50.2 
33 

701 ± 36 



1.5 



129.5 ±4.4 ±2.0 
125 — i e, 

117 - 1 -° 

12226 ± 22 



• Fits performed with the use of the lower/upper 
sideband only to parametrize the background. 

• Results from fits with alternative sets of signal am- 
plitude contributions that give equivalent Dalitz 
plot descriptions and similar sums of fractions. 

• Fits performed on a sample of 100, 000 events se- 
lected by applying a looser likelihood-ratio crite- 
rion but selecting a narrower (±la D +) signal re- 
gion. For this sample the purity is roughly the same 
as for the nominal sample (~ 94.9%). 



B. Comparison between Dalitz plot analyses of 
D+ K+K--K+ 

Table IVIII shows a comparison of the Dalitz plot fit 
fractions, shown in Table IVl with the results of the anal- 
yses performed by the E687 [8j and CLEO [9j collabora- 
tions. The E687 model is improved by adding a /o(1370) 
amplitude and leaving the i^*(892)° parameters free in 
the fit. We find that the K*{892)° width (Eq. 130]) is 
about 3 MeV lower than that in Ref. [10]. This result is 
consistent with the width measured by CLEO-c collabo- 
ration (r^, (892)0 = 45.7 ± 1.1 MeV). 

What is new in this analysis is the parametrization of 
the K+K- S-wavc at the K+R- threshold. While E687 
and CLEO-c used a coupled channel BW (Flatte) ampli- 
tude [131 to parametrize the / (980) resonance, we use 
the model independent parametrization described in Sec- 
tion lVBl This approach overcomes the uncertainties that 
affect the coupling constants and qkk of the /o(980), 
and any argument about the presence of an a (980) me- 
son decaying to K + K~ . The model, described in this 
paper, returns a more accurate description of the event 
distribution on the Dalitz plot (x 2 /v = 1.2) and smaller 
/o(980) and total fit fractions respect to the CLEO-c re- 
sult. In addition the goodness of fit in this analysis is 
tested by an adaptive binning x 2 , a tool more suitable 
when most of the events are gathered in a limited region 
of the Dalitz plot. 



Finally we observe that the phase of the </>(1020) am- 
plitude (166° ± 1° ± 2°) is consistent with the E687 result 
(178° ± 20° ± 24°) but is roughly shifted by 180° respect 
to the CLEO-c result (-8° ± 4° ± 4°). 



VIII. SINGLY-CABIBBO-SUPPRESSED 
D+ K + K~K+, AND 
DOUBLY-CABIBBO-SUPPRESSED 
D+ K+K+n- DECAY 

In this section we measure the branching ratio of the 
SCS decay channel © and of the DCS decay channel 
(|3]) with respect to the CF decay channel ([1]). The two 
channels are reconstructed using the method described 
in Sec. IIIII with some differences related to the particle 
identification of the daughters. For channel © we 
require the identification of three charged kaons while for 
channel ([3]) we require the identification of one pion and 
two kaons having the same charge. We use both the D* + 
identification and the likelihood-ratio to enhance signal 
with respect to background as described in Sec. IIIII 

The ratios of branching fractions are computed as: 

B(D+ -» K+K-K+) _ N Dt ^ K+K - K+ ^ e D +^ K+K ^ + 
B(D+ -> K+K-TT+) N D +^ K+K - n+ e D +^ K+K . K+ 

(32) 

and 

B(Df ^ K+K+tt-) ^ N d +^ k+k+ ^ e D +^ K+K - n+ 

B(Df -+ K+K-7T+) N D +^ K+K ^ + X e D +^ K+K+ ^ ' 

(33) 

Here the N values represent the number of signal events 
for each channel, and the e values indicate the corre- 
sponding detection efficiencies. 

To compute these efficiencies, we generate signal MC 
samples having uniform distributions across the Dalitz 
plots. These MC events are reconstructed as for data 
events, and the same particle-identification criteria arc 
applied. Each track is weighted by the data-MC discrep- 
ancy in particle identification efficiency obtained inde- 
pendently from high statistics control samples. A sys- 
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FIG. 12: (a) K + K~ K + mass spectrum showing a signal. The curve is the result of the fit described in the text, (b) 
Symmetrized Dalitz plot, (c) K + K~ mass spectrum (two combinations per event), and (d) the (Vj ) moment. The insert in 
(c) shows an expanded view of the 0(1020) region. The Dalitz plot and its projection are background subtracted and efficiency 
corrected. The curve results from the fit described in the text. 



tematic uncertainty is assigned to the use of this weight. 
The generated and reconstructed Dalitz plots are divided 
into 50 x 50 cells and the Dalitz plot efficiency is ob- 
tained as the ratio of reconstructed to generated con- 
tent of each cell. In this way the efficiency for each 
event depends on its location on the Dalitz plot. By 
varying the likelihood-ratio criterion, the sensitivity S of 
Df — > K + K~K + is maximized. The sensitivity is de- 
fined as S = N s I \/N s + Nb, where s and b indicate sig- 
nal and background. To reduce systematic uncertainties, 
we then apply the same likelihood-ratio criterion to the 
£>+ — > K + K~ir + decay. We then repeat this procedure 
to find an independently optimized selection criterion for 
the D+ -> K+K+it- to D+ -> K+K--K+ ratio. 

The branching ratio measurements are validated using 
a fully inclusive e + e~ — > cc MC simulation incorporat- 
ing all known charmed meson decay modes. The MC 
events are subjected to the same reconstruction, event 
selection, and analysis procedures as for the data. The 
results are found to be consistent, within statistical un- 
certainty, with the branching fraction values used in the 
MC generation. 



A. Study of D+ ->• K+ K K+ 

The resulting K + K~ K + mass spectrum is shown in 
Fig. [T27 a) . The Df yield is obtained by fitting the mass 
spectrum using a Gaussian function for the signal, and a 
linear function for the background. The resulting yield 
is reported in Table |U 

The systematic uncertainties are summarized in Ta- 
ble |Vni] and are evaluated as follows: 

• The effect of MC statistics is evaluated by random- 
izing each efficiency cell on the Dalitz plot accord- 
ing to its statistical uncertainty. 

• The selection made on the DZ + candidate Am is 



• For particle identification we make use of high 
statistics control samples to assign 1% uncertainty 
to each kaon and 0.5% to each pion. 

• The effect of the likelihood-ratio criterion is stud- 
ied by measuring the branching ratio for different 
choices. 



TABLE VIII: Summary of systematic uncertainties on the 
measurement of the Df — > K + K~ K + branching ratio. 



Uncertainty 


B(D + ^K+ K ~ K+) 


MC statistics 


2.6 % 


Am 


0.3 % 


Likelihood-ratio 


3.5 % 


PID 


1.5 % 


Total 


4.6 % 



We measure the following branching ratio: 
B(D+ -> K+K-K+) 



B(Ds 



= (4.0 ±0. ± 0.2 syst ) x 10" 



varied to ±2.5cr 



D 



and ±1.5cr 



Dp 



(34) 

A Dalitz plot analysis in the presence of a high 
level of background is difficult, therefore we can only 
extract empirically some information on the decay. 
Since there are two identical kaons into the final 
state, the Dalitz plot is symmetrized by plotting two 
combinations per event {[m 2 {K~Kt), m 2 (K-K^)} and 
[m 2 (K~K}),m 2 (K~K^)]). The symmetrized Dalitz 
plot in the D+ -> K + K~K + signal region, corrected 
for efficiency and background-subtracted, is shown in 
Fig. H"2Tb). It shows two bands due to the 0(1020) and 
no other structure, indicating a large contribution via 
D+ -> <p(1020)K + . To test the possible presence of 
Jo (980), we plot, in Fig. Il~27d). the distribution of the 
(Y]°) moment; 8 is the angle between the K + direction 
in the K + K~ rest frame and the prior direction of the 
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FIG. 13: (a) K + K + ix~ mass spectrum showing a Df signal, (b) Symmetrized Dalitz plot for — > K + K + 7v~ decay, (c) 
K + tv~ mass distribution (two combinations per event). The Dalitz plot and its projection are background subtracted and 
efficiency corrected. The curves result from the fits described in the text. 



K + K~ system in the rest frame. We observe the 
mass dependence characteristic of interference between 
S- and "P-wave amplitudes, and conclude that there is a 
contribution from — > fo(980)K + decay, although its 
branching fraction cannot be determined in the present 
analysis. 

An estimate of the 0(1020) AT + fraction can be obtained 
from a fit to the K + K~ mass distribution (Fig. 112( 0)). 
The mass spectrum is fitted using a relativistic BW for 
the 0(1020) signal, and a second order polynomial for the 
background. We obtain: 



B(D+ -» </>A H 
B(D+ - 



K+K- 
K+K-K+) 



0.41±0.08 st at±0.03 syst .(35) 



The systematic uncertainty includes the contribution 
due to Am and the likelihood-ratio criteria, the fit model, 
and the background parametrization. 

B. Study of D+ -»■ K+K+tt- 

Figure ITST a) shows the K + K + ir~ mass spectrum. A 
fit with a Gaussian signal function and a linear back- 
ground function gives the yield presented in Table [I] 
To minimize systematic uncertainty, we apply the same 
likelihood-ratio criteria to the K + K + tt~ and K + K~tt + 
final states, and correct for the efficiency evaluated on 
the Dalitz plot. The branching ratio which results is: 



B(D+ -> K+K^ 
B(Df -> K+K- 



= (2.3 ±0.3 

slat 

± 0.2 syst ) x 10" 



(36) 

This value is in good agreement with the Belle measure- 

B(Dt~¥K + K + iT' 



ment: 



B(D] 



¥K+K- 



(2.29±0.28±0.12) x 10" 3 [11 1 



Table ITXl lists the results of the systematic studies per- 
formed for this measurement; these are similar to those 



used in Sec. IVIII Al The particle identification system- 
atic is not taken in account because the final states differ 
only in the charge assignments of the daughter tracks. 



TABLE IX: Summary of systematic uncertainties in the mea- 
surement of the Df — > K + K + ir~ relative branching fraction. 



Uncertainty 


B(D + ^K+ K + tt-) 
B(Df^K+K-iv+) 


MC statistics 


0.04 % 


Am 


4.7 % 


Likelihood-ratio 


6.0 % 


Total 


7.7 % 



The symmetrized Dalitz plot for the signal region, cor- 
rected for efficiency and background-subtracted, is shown 
in Fig. ITBT b). We observe the presence of a significant 
K* (892)° signal, which is more evident in the K + ir~ 
mass distribution, shown in Fig. [TBT c). Fitting this dis- 
tribution using a relativistic P-wave BW signal function 
and a threshold function, we obtain the following fraction 
for this contribution. 

B(D+ -> K*(892)°K+) ■ B(K*(892) a -> A+tt') _ 
B(D+ -> K+K+TT-) 

0.47±0.22 stat ±0.15 syst .(37) 

Systematic uncertainty contributions include those from 
Am and the likelihood-ratio criteria, the fitting model, 
and the background parametrization. 

The symmetrized Dalitz plot shows also an excess of 
events at low K + K + mass, which may be due to a Bose- 
Einstein correlation effect [26| . We remark, however, 
that this effect is not visible in Df — > K + K~K + de- 
cay (Fig. [Hb)). 
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IX. CONCLUSIONS 

In this paper we perform a high statistics Dalitz plot 
analysis of Df —> K + K~ir + , and extract amplitudes 
and phases for each resonance contributing to this de- 
cay mode. We also make a new measurement of the V- 
wave/5-wave ratio in the </>(1020) region. The K + K~ 
5-wave is extracted in a quasi-model-independent way, 
and complements the 7r + 7r~ 5-wave measured by this 
experiment in a previous publication Q. Both measure- 
ments can be used to obtain new information on the 
properties of the /o(980) state [27t |. We also measure 
the relative and partial branching fractions for the SCS 
D+ ->■ K+K-K+ and DCS D+ -> K+K+tt- decays 
with high precision. 
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